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Abstract 

Fano profiles, which with their asymmetric character, have many potential applications in tech- 
nology. The design of Fano profiles into optical systems may create new non-linear and switchable 
, metamaterials, high-quality optical waveguides, ultrasensitive media for chemical or biosensing etc. 

CL(< In this paper, we consider an external field driven double Fano model, in which, instead of one au- 

toionizing state, there are two discrete states embedded in one continuum. We assume further that 
C , the external electromagnetic field can be decomposed into two parts: a deterministic (or coherent) 

part and a randomly fluctuating chaotic component, which is a (5-correlated, Gaussian, Markov and 
stationary process (white noise). This assumption corresponds to the case of the real multimode 
laser operating without any correlation between the modes. We solve a set of coupled stochastic 
integro-differential equations involving the Fano model with two discrete levels. We derive an exact 
formula by determining the photoelectron spectrum and compare it with the results obtained earlier 
in our previous papers. 
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. Over the last three decades, we have noticed a particular interest in the research on different ionization 

' processes of atoms in laser fields. One of the most interesting examples is the so-called laser-induced 

autoionization, on which several papers have already been published, including the study of the detailed 
characteristics of electron and photon spectra associated with photoexcited autoionization. The Fano 
model [1 for a number of discrete states and one continuum that can be diagonalized is the most com- 
monly used atomic model. Fano diagonalization, based on the Coulomb mixing of ionizing states with the 
^ , continuum, leads to a nontrivial structure of the latter [1-5]. Systems comprising autoionizing levels can 

behave in a diverse and nontrivial way, leading to interesting physical phenomena, such as quantum in- 
terferences discussed in [1-8] (and the references cited therein), electromagnetically induced transparency 
used in light slowdown 9 , and the quantum anti-Zeno effect |10| . The main feature of the Fano profile 
is its asymmetric character. Since its discovery, the Fano profile has been found in several functioned 
materials as plasmonic nanoparticles, photonic crystals and electromagnetic metamaterials. The unique 
properties associated with this asymmetric line shape lead to potential applications in a wide range of 
technologies [11]. The chemical and biological sensors are the most straightforward examples of this fact. 
An interesting review on Fano profiles in nanoscale structures is given in [12 . 

The advantage of a model comprising a Fano-diagonalized atom -I- laser light is that it can be solved 
exactly, even if relaxation mechanisms are included in the standard way. Because of the very compli- 
cated microscopic nature of the relevant relaxation mechanisms, they are usually modeled by a classical 
time-dependent random process. After this process, dynamical equations involved become stochastic 
differential equations. Except in certain special cases, e.g. chaotic white noise |13) or pre-Gaussian noise 
[14] ■ these equations cannot be solved in finite terms. 

In [3] , we treated the laser field as white noise interacting with the Fano autoionization model. Then, 
a set of coupled stochastic integro-differential equations was solved exactly. Applying the obtained solu- 
tions, we determined the exact photoelectron spectrum. In this work, we extend the formalism introduced 
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in [3] to the case of a double Fano system ^ , in which, instead of one autoionizing state, there are two 
discrete states embedded in one continuum. As it has been emphasized in [3] , the model of white noise 
for the field is interesting by itself because it describes the electric field amplitude of the multimode laser 
operating without correlation between the modes. This white noise model has been also used in our other 
paper |15| , where in a A-like configuration of electromagnetically induced transparency, the laser control 
field has the white noise component. Then the considered model is more realistic experimentally, because 
the amplitudes of the real laser light always contain some fluctuation component. As in [3 , we deter- 
mine in this paper the exact photoelectron spectrum and compare it with the results obtained earlier [3J|3]. 

The paper is organized as follows: in the second section, we describe certain aspects of our model and 
derive a set of equations for atomic operators involved in the problem. These equations are more accurate 
than those introduced in In the third section, our results are presented and discussed. Instead of 
presenting a complicated final formula, we are concerned with two interesting physical limits, i.e. for 
an infinite and finite value of the asymmetry parameter introduced by Fano [1 and generalized for the 
double Fano system in [2]. The final section presents our conclusions. 



2 A model for a double Fano system 

We start from the model shown in figure 1, where besides the ground state |0), there are two discrete 
states |1) and |2) lying above the ionisation threshold and interacting with the continuum by means of 
configurational Coulomb interaction. The two levels, as well as the continuum, are coupled to the ground 
state by a strong laser beam of frequency ojl- 




Figure 1: Simplified atomic level diagram. Configuration-interaction coupling of levels |1) and |2) to jw) 
leads to the double Fano continuum This continuum is coupled to discrete state |0) by a laser of 
frequency wl 

Following [2] we can obtain the double Fano profile describing the density of the continuum states. 
The state of this new structured continuum is denoted by round brackets |a;) (see the diagram on the 
right side of figure 1). The matrix element ri(a;) that describes the coupling of the ground state with this 
new state (the effective Rabi frequency) is given by 

r!(^) = ^=£=( ^ + + 7^, (1) 



where F = Fi -I- r2 is the total autoionization rate, the effective asymmetry parameter is expressed by 
Q = (giFi -|- q2T2)/T and a;± are the complex roots of the denominator of Eq.(l.l) given by 

^± = n + «^^, (2) 
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(3) 

where LO21 = W2 — wi, ui and W2 are the bare energies of the discrete atomic states. 
The complex amplitudes A± are given by the following expressions: 

where 

g2r2-(?iri + i(r2-ri) 
r((y + 1) 

The form equation([T]) of the radiative matrix element f2(a;) is a generalization of the corresponding 
formula of Rzazewski and Eberly [4, 5] to the case of two autoionizing levels, both of which are radiatively 
coupled to the ground state. In fact, f2(w) is a superposition of two Lorentzians and a flat background. 
Thus, for the case of the double Fano profile we have an additional Lorentzian as a result of the presence 
of an additional autoionising state. 

For the sake of convenience, we shall denote the new state by the ket jcj). As in [7, we start with the 
Hamiltonian, which describes the model with a bound state lying below the edge of the continuum, the 
continuum states and these states are coupled together by an external electromagnetic field: 

H = huoPo + J hwC^^duj + J n{uj)\0){uj\duj + H.C, (6) 

where Po and Cu,u> are occupation-number operators for the ground state and for the continuum states, 
respectively. The radiative interaction between those states is described by the function r2(w), which 
marks how strongly different points of the continuous spectrum are coupled to the bound state, and is 
given by the formula ([T]). By |0) we mean the bound state and |w) stands for the excited state in the 
dressed continuum. 

Moreover, we define the following operators 

Sc. = |0)(a.|, 

C^^, ^\u:){uj'\. ^' 

In consequence, the Heisenberg equations of motion for the atomic operators Pq, B^, Cojuj can be 
derived. They form a complete set and can be easily found by simple commutations of those operators 
with Hamiltonian (|6]). For the sake of simplicity we assume that loq = 0. 

In the Heisenberg picture, one obtains linear equations for the dynamical variables, so the equations 
for the corresponding averaged quantities can be easily found using well-known results, based on the 
theory of multiplicative stochastic processes. We assume now that has the form: 

niu:) = fiu:)iEo + E{t))e''^-\ (8) 

where E{t) is characterized by a Gaussian, Markov and stationary process (white noise), 

{{Eit)E*im^a6it-t'), (9) 

and Eo is a deterministic coherent component of the laser field. The double brackets in (jU indicate an 
average over the ensemble of realisations of the process E{t). 

Thus, we consider the following stochastic differential equation 

^ = {A + x{t)B + X* it)C}Q, (10) 
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where Q is a vector function of time and A, B and C are constant matrices. As it is known from the 
multiphcative stochastic process theory, the function {{Q)) exactly satisfies the nonstochastic equation: 

^M)l = [A + a{B,C}/2]{{Q)), (11) 

where {B, C} is the anticommutator of B and C. The result contained in equation (fTTj) is well-known in 
the literature and has been used for example in [3]. 

For the case discussed here, before performing the averaging procedure, we transform dynamical 
variables to the rotating frame: 



B^^ riu:)D^e~'''^-\ 

The new quantities D^, E^^i together with Pq satisfy the following closed set of equations: 

^ = -'1 du;\fiu;)niEo + Eit))D^ - (E* + E*it))D+], 
dD, 



(12) 



dt 
dE^ 



= liiOL - ^)D^ - tiE* + E*{t))Po + t J dio'\f{uj')\^{E* + E*{t))E^,^, (13) 
'-^i{u:- uj')E^^, + i{Eo + E{t))D^, - {{E* + E*{t))D+. 



dt 

Next, using equation pT|) we obtain a system of equations for the stochastic averages of the variables 
(double brackets have been dropped for convenience): 

^ = -aFPo~ib I dLo'\f{Lo')\\D^,-D+) + a | | dL.W|/(c.')n/(^"))l'^-'-", 
dD^ 



dt 



dt 



-ihPo+[^{uJL-UJ)-^]D^-^ J doj'\f{io')fD^,+ib I du'lfiuj'jfE^,^, (14) 
aPo + ib{D^, - D+) + tiiu - Lj')E^^, ~ ^ J dcj"|/(L^")|2(£;„.,^, +£;„„,.). 



where F = J (iw|/(w)p, b — |-Eo|- The equation corresponding to the adjoint operator £)+ is easily found 
from that for D^^ by the calculation of its complex conjugation. As we compare the formulas ([S^) and ^jp) 
in [3] with our eqns. and ((Tib ), the later contain additional terms corresponding to the coherent 

part of the laser field b which were omitted in [5j. 

Subsequently, taking the Laplace transform of equation ([Ti)) and denoting |/(a;)p — S^uj), we obtain 

zPo - 1 = -aFPo -ibj duj' S{uj'){D^, - D+) +aj J duj'du}" S{uj')S{uj")E^>u^" , 

zb^ = -ibPo + [liujL -u;)-^]b^-^ J dJS{Lo')b^, +ibj dLo'S{J)E^,^, (15) 

zE^^, - aPa + ib{b^, - b+) + i{uj - uj')E^^, duj'' S{uj"){E^'.^, + E^^^-). 

where we introduced the appropriate initial conditions for the probabilities Po|t=o = li D^\t=o = 
Euju}'\t=o — 0. As we introduce a new variable (for convenience we drop also the tilde): 

A-{z) = -ibPo -^J dJS{u')D^,, (16) 

we can write the resulting formula for D^: 

= ^ ^1 T + ^ .iJ^ ■( V / duj'S{u:')E^,^, (17) 

z + ac*/8 — — oj) z + ac* /8 — lytOL — ^) J 
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where 



B+ = 2A4 



a: 



(18) 



duj' S{uj')Ecjuj' , 



.i{T±(t)) ir±e Q-i 

and the equation for can be expressed in a similar way: 

^ A+jz) zb 

" z + ac/8 + i{u!L — ^) z + ac/8 + i{u!L — Lo) 

where A'^{z) is simply a Hermitian conjugate of (jl6p . 

To solve the whole set of our equations, we assume the separation property for Ei^i^i, i.e. 

[z - i(uj - Uj')]E^u:' = + Vio'{z)- 



(19) 



(20) 



This decomposition is of the greatest importance in the entire procedure of solving the equations in 
question. For the case of the double Fano profile, the functions (u: and rj^ satisfy the following equations, 
where we introduce oq = a/8: 



H+iz,io)Cu.^4aoPo- 
and 

H~{z,LL>)ri^ = 4aoPo + 
where 



ibA+{z) 



Z + AqC + i{u!L — t^) 

ibA-{z) 
Z + QqC* — i{ujL — w) 



4ao + 



z + GqC + i{u)L — ^) 



4an 



i?+(z,w) = 1 + aod + 
H~{z, cj) = 1 + aod* ^ 



Z + OqC* ~ i(ujL — w) 



b^d 



z + i{uj' — u) 



z — i(u)' — Lu) 



duj' (21) 



duj', (22) 



with 



4[z + aoc + i{ujL — uj)] ' 
b^d* 

4[z + oqc* — i{ujL — w)] 
S_ 1 



r V w+ — LU — iz uj- — uj — iz + 1 / 

From the analytic property of E^^r we assume that must be of the following form 



Cc.= 



ibA+{z) D+ Kcg+ 4l.+aoc+^>.-c.)] ] 



4aoPo 

H+ [z + aoC + i{ujL - ^)]H+ 



H+ 



(23) 



(24) 



(25) 



The quantity rjuj satisfies a similar equation for the complex amplitude D . By insertion the expressions 
for C^, 77^ into equations (|2ip and we obtain after simple algebra 



D 



+ _ 



4aoPo + ibgA- (z) + D' [H' {z, -iuj) - 1] 
H~ {z,—iuj) 



where 



z + aoc + — ujl) z + aoc + i{uj- — ujl) 



(26) 



(27) 



Of course, D satisfies an analogous equation after the substitution A H Using 

P^ -(| ^ and the definitions of A+ we have the following equations 



{l + aoe)A (z) 



-zbPo+'^iD^+D-), 



(1 + aoe*)A+{z) = ibPo - '°°^f ^* ( J+ + D-), 



(28) 
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where 



(jj^ — (jJl — iaQC* — iz uj^ — ujl — iaoc* — iz Q'^ + ± ^ 

T\2Lu+-iz^ 2uj^ -iz^ + 
In the similar way for Pq the first equation in (jlSp gives us 

Po ^ + ^Ke*Anz)-eA-iz)) ^1±^(,.,^ ^ ^ (30) 

z 4z 16z 

The resuhs presented above can be apphed in a steady state spectrum calculation. Such a spectrum will 
be discussed in the next section for the case of the strong external field. 



3 Photoelectron Spectrum for two Lorentzians 

As it was mentioned before, we discuss the case of the double Fano profile with two Lorentzians . For 
the model discussed here, we define the long-time photoelectron spectrum as: 

W{lu) = lim C^Ut)- (31) 

The spectrum defined in such a way can be computed directly and can be expressed in a completely 
analytical form. Thus, equations ([^5]). ([25]) and ([50]) are linear algebraic equations and their solutions 
can be easily found. The spectral distribution of excited electrons is determined from C'i^i^{t). In the 
steady state at i — > oo, only the pole at z = in the Laplace domain contributes. From the definition of 
Euiui' and its separation property, we obtain: 




Figure 2: Photoelectron spectrum for uji = lu2 = 0.5, we assume that ujl = 1-0, autoionisation widths 
Fi = r2 = 0.5 and the coherent component b = 0.1. 



cu^) = I/HP 

In consequence, the spectrum W{uj) is given by: 



Wiw) = |/(L.)n2i?eC.(0)|. 



(32) 
(33) 



Sinse, in contradistinction to the case of the single Fano profile, the analytical formula for W{u!) for the 
model discussed here is extremely involved, we do not present it here. Instead, we consider the photo- 
electron spectra for two specific cases: an infinite value of the asymmetry parameter and a finite one. As 
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Figure 3: Photoelectron spectrum for wi = W2 = 0.001, we assume that ul — 1-0, autofonisation widths 
Fi = r2 = 0.5 and the chaotic component ao — 0.12. 




Figure 4: Photoelectron spectrum for loi = 0.5, uj2 = 7.5 we assume that lol = 0.05, autoionisation 
widths Fi = F2 = 0.5 and the chaotic component ao = 0.4. 

in [2], all the frequencies (energies) are given in units of F. 

Infinite value of the asymmetry parameter (q ^ 00). For the degenerate case {u)2i = 0) the photo- 
electron spectrum is identical to that described in [3] for the case when the value of coherent component 
h is fixed, while the value of the chaotic component oq is varied (figure 2). When the values of oq are 
kept constant and the values of h are varied, the peaks are higher than in [3]. This is the case because 
equations (fMk) and Eq. ([Mb ) contain the coherent component b which was omitted from the correspond- 
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Figure 5: Photoelectron spectrum for the case of loi — 0.5, 0^2 — 7.5 we assume that ojl = 0.05, 
autoionisation widths Fi = r2 = 0.5 and the coherent component b — 0.5. 

ing equations in [3] (see figure 3). When gq is small, i.e. the coherent part of the light dominates the 
fluctuations, the photoelectron spectra exhibit a characteristic Autler-Townes splitting. 

Figures 4 and 5 present the photoelectron spectra for the nondegenerate case {0J21 7^ 0). The spec- 
trum exhibits a two-peak structure regardless of the value of b. The right peak remains almost unchanged 
while the coherent and chaotic components are varied. The left peak increases with an increase in b and 
decreases with an increase in oq. 
Finite value of the asymmetry parameter. When qi and (72 are finite, for the degenerate case (w2i = 0) 




Figure 6: Photoelectron spectrum for the degenerate case (wi = 002 = 0.5) with the chaotic component 
ao = 0, Wi = 0.5, autoionisation widths Fi = F2 = 0.5 and asymmetry parameters qi = 90, 92 = 100. 

the photoelectron spectra (figures 6 and 7) change considerably as we compare them with those discussed 
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Figure 7: Photoelectron spectrum for the degenerate case {uji = uj2 = 0.5) with the chaotic component 
flo = 0.5, lol = 0.5, autoionisation widths Fi = r2 = 0.5 and asymmetry parameters qi = 90, q2 — 100. 



W{o)) 




Figure 8: Photoelectron spectrum for oji = 2.0, uj2 = 5.0, ujl = 2.0, autoionisation widths Fi = 0.1, 
F2 = 10 and the chaotic component ao = 0.0. 

above. Figure 6 presents the photoelectron spectrum when the chaotic component is absent (oq — 0) and 
for large values of q. For the weak field case, the symmetric spectrum is revealed and it consists of only 
one peak. When the field increases its amplitude, a sharp wedge in the spectrum appears. This spectrum 
exactly reproduces that obtained by Leonski et al However, when the chaotic component is present 
(flo 7^ 0) (figure 7), the sharp wedge in the spectrum disappears - only one, symmetric peak remains in 
the spectrum and its intensity is smaller than for the case oq = 0. 

For the nondegenerate case, when the chaotic component is absent (figure 8) and for the smaller values 
of b the spectrum exhibits two peaks. However, with an increase in b, the left peak quickly decreases 
rapidly, whereas the right peak increases. As the left peak disappears, the right peak begins to split into 
two peaks. When the chaotic component is present (figure 9), the two-peak structure of the spectrum 
becomes apparent with an increase in b. This process is faster than that discussed for the case when the 
chaotic component was absent. According to the Fano theory [1], Fano zeros appearing in the photoelec- 
tron spectra manifest some quantum interferences which can occur in the system due to the fact that 
the system can achieve the continuum states via autoionization levels and by direct transitions. These 
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W(a)) 




Figure 9: The same as in Fig. 8 but for chaotic component ao = 0.5. 

interferences are weakened when the external electromagnetic field fluctuates. Therefore we can predict 
that in our model Fano zeros disappear in the photoelectron spectra. It follows from all presented figures 
that the Fano zeros really do not exist. 



4 Conclusions 

In this paper, we considered a laser-induced autoionisation model, introduced earlier in 2 , where instead 
of one autoionising state, we consider two discrete states embedded in one continuum. This model can 
be referred to as the double Fano model. Since it is interesting to simulate quantum fiuctuations by 
means of classical fiuctuations, in this paper we assumeed that the laser light is decomposed into two 
parts: a deterministic (coherent) part and white noise, i.e. a randomly fiuctuating chaotic component 
(cf 0). Subsequently, we introduced and solved exactly a set of coupled stochastic integro-differential 
equations and thus obtained the results describing the dynamics of autoionization for the double Fano 
model. The results presented here contain certain additional terms, in contrast to those discussed in [5], 
which make our findings more accurate. In particular, we derived exact analytical formulas for the photo- 
electron spectrum and compare the results with those discussed in [2, 3]. We considered the photoelectron 
spectra for an infinite value of the asymmetry parameter and, separately, for its finite value. Moreover, 
the considered spectra were discussed for the degenerate case Wi2 = and for the nondegenerate case 
(a;i2 ^ 0). We presented detailed result for both strong and weak fiuctuations regimes. As in [15], we 
believe that the model considered here is more realistic than that given in [2], because the amplitudes of 
the real laser used in experiments always fluctuate. 
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